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SUMMARY 


A  method  for  analyzing  heated  thin  wings  was  developed  under  a 
Grumman  Advanced  Development  Program  and  is  reported  in  the  paper, 

"A  Matrix  Force  Method  for  Analyzing  Heated  Wings,  Including  Large 
Deflections".^  This  paper  was  presented  at  the  symposium  "Structural 
Dynamics  of  High  Speed  Flight",  sponsored  by  the  Aerospace  Industries 
Association  and  The  Office  of  Naval  Research  in  Los  Angeles,  California, 
April  24-26,  1961. 

The  present  report  contains  some  details  on  certain  aspects  of  the 
method,  which  were  of  necessity  omitted  from  the  paper.  Several  examples 
illustrating  the  application  of  the  method  are  contained  in  the  paper 
and  might  be  referred  to  while  this  report  is  being  read. 

This  report  and  Refs.  2,  3  and  4  document  the  Grumman  matrix  force 
method  of  redundant  structure  analysis  as  developed  to  date.  The  same 
notation  is  used  throughout  these  references 0 
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SYMBOLS 


A,  P#  k,  q,  S 
Z,  Z» 


D,  F 


L,  L* 


defined  in  Ref.  (2) 

lateral  loads  induced  in  deflected  capstrips  and 
shear  panels  respectively 

length  of  a  spar  or  rib  segaent  between  panel  points 
kick  load 

constants  arising  in  systematic  calculation  of  kick 
loads 

matrices  arising  in  calculation  of  heated  structure 
flexibilities 

buckling  load  parameter 

geometry  dependent  constants  relating  the  displacements 
at  the  member  loads  to  strains  associated  with  the 
member  loads 


Subscripts* 

i*  5,  m,  n 

0,  1,  2 

or 

A,  B,  C,  D 


defined  in  Ref.  (2) 

particular  designations  corresponding  to  the  general 
subscripts  1»  J#  n  (used  in  the  derivations  to 
identify  particular  member  loads  or  panel  points) 

pertaining  to  the  Jth  member  load  of  the  structure 


pertaining  to  the  a  member  of  the  structure 
pertaining  to  the  thermal  condition 
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THEOHT 


The  basic  formulation  for  the  Grumman  matrix  force  method  of 
redundant  structure  analysis  was  developed  in  Ref*  2.  The  method  was 
expanded  later  in  company  reports  to  include  some  effects  that  are 
peculiar  to  certain  types  of  structures  or  applications.  Specif ically. 
Ref*  3  outlines  the  procedure  for  calculation  of  stresses  and  deflections 
resulting  from  non-uniform  thermal  strain,  while  Ref.  4  explains  the 
inclusion  of  Poisson* s  ratio  and  sweep  coupling  effects* 

A  recent  paper.  Ref,  1,  further  extended  the  method*  as  applied  to 
the  particular  case  of  a  thin  wing,  to  include  (1)  the  Influence  of 
in-plane  thermal  stresses  on  the  flexibility  of  the  wing  for  lateral 
load,  (2)  the  limiting  case  of  increased  flexibility  from  this  cause, 
namely  thermal  buckling,  and  (3)  the  stiffening  effect  that  results 
from  large  deflections*  The  emphasis  in  the  paper  was  on  the  basic 
formulation  and  examples,  and  due  to  space  limitations,  only  the  most 
important  formulas  could  be  given  in  some  instances* 

The  following  sections  give  a  more  detailed  development  of  the 
formulas • 


It  will  be  assumed  in  the  discussion  that  follows  that  the 
wing  structure  is  symmetrical  about  a  mean  chord  plane  and  that 
it  can  be  idealized  in  the  stringer,  shear  panel  and  shear  web 
manner*  This  is  one  of  the  standard  methods  for  idealizing  a 
multi-spaj;  multi-rib  wing  and  consists  of  representing  the  wing 
by  capstrips  carrying  axial  load  only,  interconnected  by  shear 
panels  which  can  carry  tangential  edge  loads  only.  Using  this 
approach,  one  lumps  the  covers  into  equivalent  capstrips  over  the 
spars  and  ribs;  the  contributions  of  the  spars  and  ribs  themselves 
are  calculated  and  added  in.  If  appropriate,  intermediate  shear 
lag  members  are  introduced  as  well*  Because  of  the  shear  panel 
assumption,  the  axial  loads  in  the  capstrips  vary  linearly  from 
one  panel  point  to  the  next. 

Based  upon  this  idealization,  the  member  loads  in  the  statically 
determinate  structure  and  the  member  flexibilities  can  be 
determined  by  the  procedures  already  given  in  Refs*  2,  3  and  4» 
Following  these  procedures  further,  one  can  determine  the  stress 
distribution  and  the  deflections  due  to  applied  loads  and/or 
temperature  changes  under  the  assumptions  of  linear  behavior  of 
the  structure » 
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The  remaining  sections  are  concerned  with  the  procedures 
which  are  applicable  in  cases  where  in-plane  and  out-of-plane 
components  of  load  and  deflection  cannot  be  directly  superposed* 

C*  Interaction  Problem  -  Small  Deflections 

The  problem  of  analyzing  a  wing  under  an  arbitrary  temperature 
distribution,  neglecting  the  interaction  effect,  can  be  handled  by 
considering  the  wing  to  be  subjected  to  two  component  temperature 
distributions,  symmetrical  and  anti-symmetrical  about  the  mean  chord 
plane  respectively,  which  when  superposed,  add  up  to  the  given  distribution* 
The  wing  is  analyzed  under  each  separately,  and  then  the  resulting  stresses 
are  added  algebraically*  The  temperature  distribution  symmetrical  about 
the  mean  plane  Will  give  rise  to  stresses  symmetrical  about  the  mean  plane 
and  displacements  in  the  plane,  while  the  anti-synnetrical  temperature 
distribution  will  produce  anti-symmetrical  stresses  and  displacements 
normal  to  the  plane.  The  latter  add  directly  to  the  stresses  and  de¬ 
flections  caused  by  lateral  applied  loads. 

An  obvious  assumption  in  the  foregoing  procedure,  and  a  valid  one 
in  many  cases,  is  that  the  two  analyses  can  be  treated  independently. 
However,  if  the  in-plane  stresses  are  large  enough,  there  can  be 
noticeable  interaction  between  the  two  thermal  stress  distributions 
themselves  and  between  the  in— plane  stresses  and  the  lateral  displacements 
due  to  applied  loads.  In  this  case  the  flexibility  influence  coefficient 
matrix  for  the  structure  may  be  appreciably  altered. 

In  the  development  of  the  method  for  accounting  for  this  effect  it 
is  assumed  that  the  deflections  are  small  in  the  sense  that  no  membrane 
type  stresses  are  developed  due  to  deflections.  The  two  Idealized 
structures  used  in  the  symmetrical  and  anti-symmetrical  analyses  for 
a  given  case  will  be  geometrically  identical.  Element  flexibilities 
however  may  be  different,  inasmuch  as  the  stresses  in  one  case  depend 
upon  the  bending  and  torsional  flexibility  and  in  the  other  upon  the 
in-plane  flexibility*  With  the  "cold”  structure  flexibility  influence 
coefficient  matrix  and  the  symmetrical  and  anti-symmetrical  thermal 
stress  distributions  already  calculated  by  the  methods  given  in  Refs.  2 
and  3,  the  two  analyses  can  be  coupled  by  the  procedure  to  follow. 

In  the  development  of  the  coupling  formulas  it  is  convenient  to 
consider  first  the  symetrical  thermal  stress  distribution*  in  this  case 
the  loads  in  pairs  of  upper  and  lower  capstrips  and  shear  panels  are 
identical.  This  load  distribution  is  self  equilibrating  only  as  long  as 
the  structure  remains  undeflected  laterally. 
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When  the  spars  and  ribs  bend  and  the  shear  panels  become  warped#  resultant 
forces  normal  to  the  plane  are  induced*  To  illustrate,  consider  the 
segment  of  spar  or  rib  shown  schematically  in  Pig.  1.  As  shown,  the 
average  slope  between  panel  points  is  taken  as  the  element  slope.  To 
calculate  the  magnitude  of  the  lateral  force,  which  will  be  called  a 
kick  load  and  denoted  by  Z  ,  consider  vertical  components  of  forces  at 
the  panel  point.  Summation  of  these  forces  gives 


U-O 

l 


a) 


where  Zc^  t  is  the  total  load  in  the  upper  and  lower  caps  trips;  Sa  ,  §x 
and  6Z  are  the  deflections  at  three  conseetive  panel  points,  Mid 
?,  and  'lz  are  the  distances  between  them. 


Fig.  2  shows  a  pair  of  deflected  quadrilateral  shear  panels, 
together  with  their  projection  on  the  undeflected  mean  chord  plane*  The 
mean  positions  of  the  four  corners  are  identified  by  the  four  panel  point 
deflections  dA  ,  ,  6C  ,  and  60  *  The  In-plane  shear  flows  lcA6  a  , 

etc*  have  a  lateral  resultant  which  is  called  H'  and  can  be  calculated  by 
again  summing  forces  vertically.  Thus, 


Z  -2 CJ,  ( +■  kAD )  ( kAP  -v  kBC )  cffe 

4-  (  k6C  -v  kco )  <fc  -  i  k*,  +  kcp )  ZD  j 


(a) 


In  the  case  of  the  shear  panels,  the  force  in  eqn.  (2),  which  is 
assumed  to  act  at  the  intersection  of  the  diagonals,  is  replaced  by  a 
statically  equivalent  set  of  forces  at  the  four  comers  of  the  panel* 


In  the  formulas  to  follow,  the  matrix  notation  of  Ref.  2  will  be 
extended  and  used*  The  kick  loads  will  be  written  in  column  matrix 
arrangement  and  designated  {P^}.  They  are  dependent  upon  the  deflections, 
which  in  column  order  are  called  \  <Tn} ,  and  the  in-plane  member  loads  which, 
for  convenience,  are  written  in  diagonal  matrix  form.  The  latter  are 
designated  The  formula  from  which  the  kick  loads  are  obtained  is 

then 


(3) 
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where 


F  =  P  -a  r 

*  rT r,  V-Ai.i  vi3  '^•1 


The  and  ICKA  matrices  are  made  up  of  geometry  dependent  constants 
derived  from  eqns.  (l)  and  (2).  They  can  be  obtained  by  a  systematic 
procedure  which  is  discussed  later* 

Recalling  from  Ref.  3  that  the  deflections  due  to  lateral  loads  and 
anti-symmetrical  thermal  strains  are 


Aj  jPnl  +  \  ii 


the  total  deflections  may  be  written 


..  _  Amij  ^  |  1  Ame,  '  Awr.  Fj 


Solving  for  l  Snl, 


ll-WWW +  [af  W 


(5)  • 


where 


[ClJ=(:I]  - .  W  [Fvlg-.  ■ 

It  is  seen  that  tp^3  '  LAmn  J  is  the  required  matrix  of  flexibility 
influence  coefficients.  Multiplication  of  [.Aj  bylD^y'introdueea  the  ehange  in 
flexibility  due  to  in-plane  thermal  stresses,  while  multiplication  of 
hylDJ"' takes  care  of  the  interaction  between  the  symmetrical  thermal 
stresses  and  the  thermal  displacements  discussed  previously. 

With  the \ 5 n\ known,  it( is  now  possible  to  substitute  in  eqn.  (3)  to 
determine  the  kick  loads  {P^  With  these  known,  the  total  anti-symmetrical 
components  of  the  member  loads  may  now  be  determined  as 
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In  eqn.  (5),  the  subscript  "as*  refers  to  the  lateral,  or  anti-symmetrical , 
analysis*  The  +  and  -  signs  indicate  the  top  and  bottom  surface  member 
loads  respectively*  In  those  cases  where  the  section  properties  are  the 
same  for  the  anti-symmetrical  and  the  symmetrical  analyses,  the  in-plane 
loads,  \ may  be  added  directly  to  the  loads  given  by  eqn.  (7)*  If 
this  is  not  the  case,  then  it  is  necessary  to  calculate  actual  stresses 
for  the  two  cases  and  then  superpose* 


D.  Analysis  -  Systematic  Calculation  of  Jk.  and  Cj^ 

In  order  that  the  method  be  appropriate  for  routine  use,  the  entire 
set  of  kick  loads  {P*j^  must  be  obtainable  by  a  systematic  procedure  for  any 
given  in-plane  thermal  stress  distribution*  The  result  of  one  such  method 
is  indicated  in  eqn*  (4)* 

Consider  first  a  typical  segment  of  spar  or  rib,  as  illustrated  in 
Pig.  (3).  In-plane  member  loads  q.  and  are  assumed  to  act  on  this 
particular  pair  of  caps  trips;  deflections#-  and  <!L,  together  with  the 
corresponding  kick  loads  PI  and  PI  are  shown  as  Well*  Referring  now  to 
eqn.  (1),  the  contributions  of  this  pair  of  capstrips  to  PJ  and  P^  are 


4 


"A 


t 


p:  = 

0  n 


<fc  ~ 

L  " 


These  expressions  can  be  written  in  matrix  form; 


R 


*  . 


p: 


2 


/it  -  A 

k] 

m  , 

The  contributions  of  this  pair  of  capstrips  to  the  B,  and  C,  matrices 
are  thus;  mi  in 


[§J  = 


1 

A  -w 

-;>6 . - 

r 

y* 

-  . . ^ . . 

p" 

%  "" 

C 

Q  - 

A 

-y\ 

§ 

u 

Ref  erring  now  to  the  typical  cover  shear  panel  shown  in  Pig*  (4),  the 
in-plane  member  load  q. ,  in  combination  with  the  lateral  deflections  of  the 
four  corners  6 ^  dL, at  and lances  a  kick  load  Z»,  a*  described  by 
eqn.  (2)*  Distributing  half  of  this  load  to  each  pair  of  diagonally 
opposite  comers, 

p*  ■  — -  rj  —  M 

.  *  2(viU  K  ~  w. TiU"; 

p*  —  4  . ZZ1 


Mm 


41 
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Substituting  for  Z*  from  aqn.  (2)  and  rewriting  in  matrix  form. 


p; 

p; 
p: 
p' 

-  o 

The  con' 


Mu 
U* 
Mi. 


V 

bribu^lons  of"  this  pair  of  cover  panels  to  the  B  .  and  C. 
matrices  are  thus?  111 


>? 


t  stc) 
*  L) 


L 


4 

A 


[B  .] 

w  ana J 


i 

A 

6 

C  '• 

0 

j- 

lc,]  = 

t MCI 


A _ d 

& 

c  • 

D 

i 

.  X 

The  B  .  and  C_.  matrices  can  easily  be  formed  in  this  manner  for 
every  pair oT  capstrips  and  for  every  pair  of  cover  shear  panels  in  the 
structure*  They  may  then  be  added  together  in  the  same  way  that  member 
flexibilities  are  added  to  form  the  member  flexibility  matrix 
of  Ref.  J,  to  yield fc^Tand  tC -lol  • 

It  is  recommended  that  in  numbering  the  member  loads  the  numerical 
subscripts  designating  member  loads  acting  upon  any  one  member  be 
reserved  uniquely  for  that  member*  (For  example,  even  though  axial  loads 
in  adjacent  segments  of  capstrips  must  necessarily  be  equal  where  they  join 
together,  it  is  recommended  that  the  member  load  numerical  subscripts  be 
selected  such  that  they  do  not  repeat.)  This  is  desirable  for  several 
reasons*  In  this  instance,  it  is  beneficial  because  vixen  the  B  ,  and  C  . 
matrices  for  all  of  the  members  are  collected  into  the  B  .  and  matrices, 
there  will  be  no  overlapping  of  individual  elements*  Tfcal  in  turfi  means 
less  chance  for  error* 


E*  Thermal  Buckling 

It  is  well  known  that  the  increased  flexibility  of  structures  for 
lateral  load  is  related  to  the  proximity  of  the  in-plane  load  to  buckling. 
The  same  matrices  used  to  calculate  the  heated  structure  flexibilities  can 
be  used  to  determine  the  factor  by  which  a  given  in-plane  loading 
distribution  must  be  increased  to  cause  buckling. 
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If  the  structure  Is  subjected  to  in-plane  loads  only ,  then 
according  to  eqn.  (5)* 


R}  =  [a„]  [fJ  \i 


Examination  of  eqn.  (8)  indicates  that  equilibrium  configurations 
other  than  the  undeflected  case  may  be  possible*  If  they  are,  then  the 
deflections  cfn  and  the  loads  entering  CFmnl must  satisfy  eqn.  ($)  and 
can  be  determined  from  it*  Assuming  that  the  in-plane  loads  increase 
proportionally,  then  a  given  load  distribution  becomes  at  buckling 
aC'ai83*  a  is  the  factor  which  when  multiplied  by  each  of  the  ins-plane 
loads  gives  the  load  distribution  causing  budding*  Eqn.  (6 )  can  now  be 


written  as 


s\  =  [a  1  If  V  \s 

^nj  [’  mnj  [  un 


which  is  recognised  as  an  eigenvalue  equation*  The  lowest  value  q 
satisfying  this  equation  defines  the  in-plane  loading  condition  for 
buckling.  The  corresponding  eigenvector  1  <f)A}gives  the  buckled  shape. 


An  implied  assumption  in  linear  wing  analysis  methods  is  that  the 
lateral  deflections  are  small  enough  so  that  the  attendant  straining 
of  the  middle  plane  is  negligible.  It  is  quite  possible  that  this 
assumption  may  not  be  valid  for  the  case  of  a  thin,  low  aspect  ratio 
wing.  Chi  the  contrary,  as  a  consequence  of  large  deflections,  there  may 
be  a  stiffening  of  the  wing  due  to  the  development  of  in-plane  membrane 
type  stresses  which  help  carry  the  load* 

In  Ref.  l  a  method  is  developed  for  including  this  effect.  It 
permits  the  determination  of  the  deflected  shape  for  a  given  set  of  applied 
loads,  together  with  the  attendant  stress  distribution,  and  the  flexibility 
of  the  structure  for  small  displacements  about  the  deflected  shape.  The 
latter  is  stated  by  means  of  a  tangent  flexibility  matrix  which,  of  course, 
is  dependent  upon  the  particular  loading  condition* 

The  first  step  in  the  procedure  is  to  calculate  the  stress 
distribution  in  the  middle  plane  consistent  with  an  asstned  or  calculated 
approximate  set  of  deflections.  This  calculation  is  similar  to  the 
calculation  of  in-plane  stresses  due  to  non-uniform  thermal  expansion!  in 
the  present  case  the  strains  are  due  instead  to  deflection. 
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Formulas  for  the  strains  in  the  statically  determinate,  stress  free, 
deflected  middle  plane  are  given  in  Ref «  1  for  spars  and  ribs  and  for 
rectangular  shear  panels.  To  complete  the  collection  of  formulas  that 
may  be  needed  for  swept  wing  analyses,  expressions  are  derived  in  the 
Appendix  for  shear  panels  of  general  quadrilateral  shape. 
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Ref.  3  develops  the  following  expression  for  the  virtual  work  done 
in  a  quadrilateral  shear  panel? 

p 

r5r(2e;bnX+?ss,)dVMi(Ll(5x.i:i,  7  1  <“> 

iJV  •'  ■ 

The  various  terms  of  this  expression  are  explained  in  Ref.  3 1  in  particular 
a  is  the  member  load  shear  flow,  L„  and  L\,  are  functions  of  the  shape 
of5 the  panel,  and  5  and are  the  direct  and  shear  strains  associated 
with  a  line  bisecting  the  Ingle  between  the  two  opposite  sides  of  the 
quadrilateral  that  are  most  nearly  parallel*  See  Fig.  5*  It  is  the 
primary  purpose  of  this  section  to  derive  equations  for  the  strains  ex  and 
K  when  they  arise  due  to  large  deflections. 

2*  Derivation  of  Expressions  for  e  and  Tf 


Referring  to  Fig.  5,  the  bisector  of  the  angleflf-  oOis  designated 
the  x  axis;  the  bisector  of  the  angLe($~/$}  is  called  the  n  direction. 
Then  in  terms  of  the  lateral  deflections  of  the  four  corners  6 6 6  q 
and  <5n»  the  average  slopes  ^  w  and  d  w  are? 

b*  dn 


dw  .  t 

■  T .  35 

OX  z. 


4  -4 


±  if  4  -  4  jl  4  “4 

dn  ~  a  C  *  ~nrr 


The  quantities  X  X  gn#  JGD  and  XaD  are  the  lengths  of  the  four  sides 
of  the  quadrilateral.  Use  of  eqn*  ^(a),  page  34  of  Ref*  5  now  yields 
for  the  slope  in  the  y  direction: 


Sg-=-gf  COti(»  +«  +  *♦•<$) 

CSC  4r(<x4-fl  +  ?*/•)’ 
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The  strains  ex  and  Tf  due  to  large  deflections,  whose  finite  difference 
equivalents  for  rectangular  panels  are  given  on  pages  1$  and  16  of  Ref.  i» 

become  r  A  ai 

^  L  J  i  &WY  ] 


nm 

=  ~  -37]  esc  xG 


Substituting  eqns.  (2A)  and  (3A)  in  the  preceding  gives  the  desired 
equations. 

3«  Equivalent  Strain  If 

— . — - 

* 

The  right  hand  side  of  the  virtual  work  expression,  eqn.  (1A), 
consists  of  two  terms,  one  in  e  and  the  other  inY.  Both  strains  are 
functions  of  the  same  panel  displacements  d\,  cfR*  d>»  6  ^  however,  and 
it  is  therefore  convenient  to  combine  everything  into  a  single  term.  This 
can  be  done  in  the  following  way.  An  ^equivalent*  shear  strain  f  is 
defined  such  that  eq 


c\,  Xo^  "  ^  +  ^11 


\  €x  *  X 


^  18  9een  io  be  a  quadratic  function  of  JA,  6b,6q,  <5d,  as  follows* 

Tt=-U«4441  fa„o„  4*^1  A  444  4  \ 


Where 


j-  i  ^  +iUh6) 

CSC  7U  +  +  <T) 
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